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Abstract 

Fermion determinant is computed analytically on extremely large 
lattices N T — > 00 in the toy model approximation in which action is 
truncated so that in the Hamiltonian limit of a T —* all terms of order 
a T /a a are discarded. Chiral condensate is studied in the area of small 
( m « T) quark masses. 
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1 Introduction 



Finite-temperature studies of the fermion contribution in LGT have been 
on the agenda over a number of years and both experiments and analytical 
models were tried to challenge the problem. A crucial issue for such studies 
is the fate of chiral symmetry in continuum limit. Urgent need to perform 
calculations in full QCD motivated substantial efforts made to overcome 
technical problems caused by full dynamical treatment of fermions. Recent 
developments in the technology of lattice gauge theory made the computa- 
tion of fermion contribution feasible. Calculations with dynamical fermions 
and small lattice spacing, however, are still almost unthinkably expensive. 
Therefore even crude models that allow an analytical solution might be use- 
ful. 

Unfortunately we can hardly attack the full problem, hence the approxi- 
mations which hopefully capture some of the essential features of the physics 
may be considered. Here we attempt to study such a problem in 'electric 
toy' model approximation 0, where all terms proportional to l/£ has been 
removed []. Several of the tools we employ exist in some form in the earlier 
literature on the subject. For example in pure gluodynamics this approxima- 
tion permits to discard the magnetic part of the action EL In this aspect 
such approximation is similar to strong coupling approximation. Moreover, 
effective action formally coincides with obtained by 0, P in g 2 » 1 limit. 
On highly anisotropic lattice £ >> 1, however, the condition g 2 » 1 is su- 
perfluous and for N T » 1 weak coupling region is penetrable which allows 
to study continuum limit in this model. 

In our previous article |J pure QCD was considered in such approxima- 
tion. On extremely large lattices ( N T — > oo ) the model can be solved 
analytically and allows to study running coupling behavior in the continuum 
limit. Results obtained for the Callan-Symanzik beta function are not too 
surprising, but hardly close to reality since they predict trivial asymptotic 
freedom g 2 ~ a T . The simplest way to move toward a more realistic case 
of interacting quarks and gluons is to study dynamical fermions in the same 
pattern. Although it scarcely makes the model entirely realistic , it seems 
not incurious to study the chiral condensate behavior at finite temperatures 
within mentioned approximation. 

1 The anisotropy parameter £ is defined as £ ~ a a /a T , where a T and a a are temporal 
and spatial lattice spacings, respectively. 
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In this paper we focus our attention on lattice formulation of QCD using 
Wilson fermions, only for convenience reasons. As it is known both Wil- 
son and Kogut-Susskind formulations have their advantages (comprehensive 
analysis of pros and contras is given in |||7|]) and, of course,the results should 
agree in continuum limit. Due to well known "no-go" theorem || there is no 
straightforward way to remove fermion spectrum degeneracy without break- 
ing chiral invariance. For Wilson fermions chiral symmetry is broken ex- 
plicitly by the Wilson term needed to remove the unwanted doublers on the 
lattice. This makes it difficult to study the spontaneous breaking of chiral 
symmetry using Wilson fermions. Staggered fermions, indeed, provide mass- 
less fermions on the lattice but the continuum notion of chirality is not well 
defined for staggered fermions due to their inherent lattice construction M. 

For the fermionic part of action we choose the form suggested in Jl(| (see 
also 0) 

\ x n=l x / 

where n/ is the number of flavors^] and 

D x'x = 2 (x) S x ,x>-n + T y^ ^l {%') ^x,x'+n ~ r5 x , )X (2) 

with 

D x'x = 2 7 ^ o {x) S x ,x'-o + T + 2 y ° U o ( x> ) $x,x'+o - (ma T + r) 5 X > >X (3) 

Extra term with finite Wilson parameter r (0 < r < 1 to guarantee positivity) 
becomes an irrelevant operator for ordinary fermions in the infrared limit, 
while being a relevant operator for mirror fermions in the high-energy regime, 
in fact , generating an effective mass for mirror fermions^ 

Inasmuch as in suggested approximation we omit the terms of l/£ order, 
the 'toy' action will be simply a 3 X^z-^L'Vv- Since 'toy' fermions propa- 
gation in space direction is suppressed by factor l/£, they may hardly claim 

2 Later (unless expressly specified otherwise) we shall consider rif = 1, because the 
generalization is evident. 

3 As it is shown below, in a given approximation it is equal to m e f / = tn + ^- arctanh r 
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to be dynamic in the proper sense of the word, nevertheless they bring a 
nontrivial term into the effective action, which we try to calculate. 

Fermion fields ip x can be presented as a linear combination of creation 
at (annihilation a PjCr ) operators of the particles and antiparticles h^_ p _ a 
(b_ Pi _o.) with four momenta p and projection of spin a 

p,a 

1p x = ^2 {Up**v,o + U-p-ob-p,-<r) ■ (4) 



p,cr 



In a standard quantization procedure ( see, e.g. [f[^] ) such operators serve as 
integration variables, however, their combinations may be used as well, on 
condition that corresponding Jacobian does not turn to zero. To choose a 
convenient combinations we remind that in a standard representation 

70= (J -l) ; 7 » = (-a n ^") ; 75= (-l "o 1 ) (5) 
projectors 

l + 7o /l 0\ l-7o /0 fl\ , . 

— = o o J ; — = V 1 J (6) 

divide bispinors ip into two components ip^, each including only one two- 
component spinor 

/(+) l±7o, fl> W \ ,(-) 1-70/ / \ 

(similar to free particles and antiparticles at rest) and are the eigenvalue of 
inversion operator 

P^(±) = ±i^(±) (8) 
Taking into account that (5 z _/j )a; /) = b x ,x-{i' = &x+p,,x'i an d presenting 

r±7o _ l + 7 r±l r =F 1 1 - 7o ^ 
2 2 2 2 2 U 

we may rewrite (Jj) in terms as 

- ^ = VT^ (^f^A +) + ^A^-A ; (10) 
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where 



K x > = <$xx' I ^-^^f ,t-i ° - m'<k.t ) (11) 



with 

a T m r = arctanhr (12) 

and 

, ma T + r 



r 
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ma T cosh (a T m r ) + sinh (a T m r ) = (m + m r ) a r + O (a^) 



(13) 

By gauge transformation all £/o (x, i) matrices may be diagonalized si- 
multaneously: U {*,t) au = 8 au U {x,t) aa , therefore A xx , in ([TT]) is simply a 
set of N matrices N T x N T . If we fix the static gauge Uo (x, t) aa = uj (x) Q , all 
such matrices will depend only on £ — t' therefore they may be diagonalized 
by the discrete Fourier transformation exp (ijfptj with integer k in for pe- 
riodic border conditions on fermion fields at time direction, and half-integer 
k in for antiperiodic ones, demanded by Osterwalder-Schrader |nj positivity 
condition?]. Formally a solution may be obtain with the help of the simple 
equation 



IB -t+lN T 
Uu 



1 / ^illh. 1 i^Lk\ U -i2*ht ( n \ f 1< 

— }_Ue »», exp'** J e ^ I n _ WiVr J(-l 

r fc ^ ' «=-oo \ 2 / 

(14) 

where I? = for periodic border conditions and B = 1 for aperiodic ones. 
The summing over / is the tricky procedure and regularization is evidently 

needed. If, however, we simply change exp ^— i^- J — > exp (— i<p) and J2k " 

Jo ^ %t ' ^ mean l — in ( |I~4"D so that the difference between periodic and 
antiperiodic border conditions is lost, which is regarded as very essential in 

4 Osterwalder and Schrader jlj] developed a mathematical procedure that allows the 
reconstruction of Hamiltonian and physical Hilbert space from continuum field theory 
defined in Euclidean space. For Wilson formulation of lattice gauge theory this condition 
was shown to hold by Osterwalder and SeilerflfjJ. 
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A 



this case. Hopefully, the computation of determinant 

/ m'a T At/J(x,0) 

-XU (x, 0) m'a T \U ] (x, 1) 

— XUq (x, 1) m'a T 





^ A(-l) B E4(x,JV T -l) 



-\(-l) B U (x,N T 



At/ t (x,iV r -2) 
m'a T 



(15) 

where A = - — ^— ^ , in the considered extremely simple case can be done 
straightforwardly and, incidentally, this does not even need the gauge fixing 
(see Appendix). After the integration over ipx fields we get 



- S e J f = In det A f A = ^ In det A^A a = ^ H* 



with 
where 



[m 



, AT -L wir . X mr 

M Nt + (-lf T -e—tt a + -e- — 



(16) 
(17) 



cosh (N T arcsinh m'a T ) ~ cosh (^) + (a 2 ) ; N T = 2k; 
sinh (iV T arcsinh m'a T ) ~ sinh (^H + O (a 2 .) ; N T = 2k + 1. 



Therefore for ei>en N T we may write 

-S^ = 53 In (V + cos 2 (p a + {-if (^cosh ^ 



m + 2m r 
T 



(18) 

+ cosh (t^-)^ cos<p a 
(19) 



and , respectively, for odd N T 

N f f 

— S e J^ = 5^ m ( P 2 ~ COs2 'Pa + ( — 1) B ( C0Sn (jr) _ COSn 
a=l ^ ^ 



m + 2m r 
T 



where 



2 i 2 "I - ^^7* i 9 

p = 1 + sinh — h sinh — . 



Both ( |T9"D and ( p0|) can be nicely expressed as 



N 

S e J^ = In J^J (z r — COS<p Q ) ^2 — (— l)^ COS<p Q j 



cos<p c 
(20) 

(21) 



(22) 



withQ 

z = (-l) cosh—; z r = (-1) cosh (23) 

In case of SU (2) gauge group (pi = — y? 2 = § = § arccos |, where x is the 
character of fundamental representation, so (B^) can be easily expressed in 
invariant form 

- S^' = in n (-h £ + (-!)«* |) (cosh + (-1)' f ) ; 

(24) 

In case of SU(3) gauge group we can express S e J* through the characters of 
fundamental representation \ = E«=i e%iPa with the help of simple relation 

3 

I [ (l + pe iip<y + qe~ iipa ) = 1 + p 3 + q 3 - 3pq + (p + q 2 - 2p 2 q) x + 

a=l 

(q + p 2 - 2pq 2 ) x* + pq 2 X 2 + P 2 QX* 2 + W 
that, e.g. in case of even N T , leads to 

- Sf = In / (cosh + In / (cosh ^±^) ; (26) 

with 

/ (x) = x 3 PQ 2 + *±^x 2 PQ + (-if X * + g (x * )2 + ^Cr; (27) 

and 

P(x) = l + t^L; Q (x) = (-if - ±- (28) 
x 2x 

As it can be seen from (|26|) the effective action S F looses Z (N) invari- 
ance, already in zero order in Quark effects lift Z (N) degeneracy and 
only the phase in which Polyakov loop is real is stable. Other phases are 

J The parameter r was introduced to shift the mass of mirror fermion. As it follows 
from definitions ( |l2| ) and ( |l3| ) m! — m r = m and ml +m r = m + 2 m r , so the suggested 
parametrization r = tanh a T m r , makes the machinery of such shift more transparent. By 
increasing m r one can make the doublers as heavy as the cutoff. It is easy to check that 
if Imr = it does not introduce any imaginary part into the action. 
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metastable. As it can be seen from P^D, the strength of Z (JV) symmetry 
breaking decreases with increasing mass 



r. e ff n ^ T f m + m r , ^\ / ->\B ( , ->\N^ m m+2m r \ Rev 
Sjp ~ 2iV f — - In 2 J + (-1) [{-I) T e~- +e ~~" * 



T 



(29) 

in agreement with the reasons given in ]15[ . This will modify the thermody- 
namic properties of metastable phases, the lifetime of which increases with 
decreasing strength of Z (N) —symmetry breaking. 

It is nameworthy to note that, if the integrand of partition function is 
highly peaked in proximity to x min = J2 a cos<^™ in , with fixed y?™ 111 , the par- 
tition function has zeroes in complex m-plane at 

m a = ±i (ip a + irn) T r ; m' a = ±i (ip a + nn) T+2m r ; n = mod (B + N T + 1) 

(30) 

As it can be seen from (|18"D , M^ T functions are polynomials^ in m'a T of 
order N T , therefore, at finite a T the equation S a (m') = has N T solutions 
{m' t } which are, generally speaking, different and formally can be written as: 

m > j JYrSinh^axccosh^^p^ ( e ^ n « + e "^?t))) ' N T = 2k; 
Y= \ iV.sinh^arcsinhtipi^e^^ + e-T^)); N T = 2k + l. 

(31) 

For finite a T and N T the zeroes of det At A dispersed in m-plane, but with 
a T — > 0,and N T —>■ oo they steadily concentrate around m a and m' a . 

From fl2"4|)it is evident that fermion loop contribution increases in the 
parameter area where such zeroes approach the origin of coordinates in m- 

l , 2 

plane (while T is fixed) and simultaneously py 2 — > 1. 

As it is known (see e.g., (JTHJ and references therein), the chiral condensate 



e expressed m terms of its zeros m t = — — v as: 

1 L cosa(m r a T ) a T 



-Y^—. (32) 
V a ^— ' m—m t 



where V a is spatial lattice volume. The behavior of chiral condensate in 



QCD can be investigated |L6| having studied the distribution of partition 



Mn t functions are very similar to Chebyshev polynomial, e.g. 
cosh (n arcsinh \J z 1 + l) = T n (x) — cos (n arccos x) 
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function zeros in the complex quark mass plane. Zeroes at (|3T)|), certainly, 
don't exhaust the whole set of {m t } , but seem likely to represent a very 
important part of it. 

To finish with this issue we want to note that in a toy model approxi- 
mation the fermion contribution is simply a 'potential'-like term and can be 
totally incorporated into the measure. For example in case of SU(2) gauge 
group the modified measure has the following form 

d\i — ^cosh — + (-1) cos— J I cosh — h (— 1) cos — I sm — dtp 

(33) 

For arbitrary flavor number n/ it can be written as 



m . i -,sB+n t ¥\ nf ( ,m + 2m r b ip 



v 



dfinj = (^cosh— + (-1) Tcos ^J [cosh h(-l) C0S 2j sin 

(34) 

Therefore, in a toy model approximation the fermion contribution does not 
seriously complicate the partition function Z and thereby the computation 
of any average values . Here we concentrate on the evaluation of = 



2 



T dlnZ 
Va- dm 



2 Chiral condensate 

In case of SU(3) gauge group in the area ^ >> 1 one may write 

/^\ ~ OIL + -^L + ( (~ 1 ) JYt 9z + / x \ + 

^ ' z dm z r dm I z 2 dm zl dm I 

( m m + 2m r 

~ tanh — + tanh — 

V T T 




T 



which means strong breakdown of chiral symmetry. Therefore, gluon envi- 
ronment influence on (V>V) value became inappreciable when s infinitely 
increases. 

To compute (ipip) in the area ^ << 1; ^ << 1> we use 'mean spin 
method' |I7j which in fact is a modification of the method of quasiaverages 



i; 



-<A^= JA5 Un*-J^cos Vx j JI^x (36) 



The effective action S = S (x) can be computed as 

X n=0 x 

and in the lowest order can be obtained simply by substitution 2 cos <^ x — > 
X- In such approximation the results coincide with a mean field method. 
Ordinary averages may be evaluated by integration over x- 

Let us consider for simplicity the case of SU(2) gauge group were the 
gluonic part of action for a 'toy' model |lj can be written asQ 

- s = J2 cos v cos ^ = p E ( 38 ) 



and fermion contribution is incorporated in the measure dp, given by fl33| ) . 

Since we introduced the collective variable x — ~m J2 X cos V?x the corre- 
sponding Jacobian -< 1 y% must be computed. Having written 

S (§Aj - £ cos = £J exp { 3 (f - £ cos ^ } £ (39) 

we find 

-< 1 ^= exp (-J\#fi) = exp {N d a (f s - L («)) } ^ (40) 

where^] 

r-{-?}*-(*-5)(-5)^ < 4i » 



/■47T 

-L(s) 



or 



2 6 \ /i (s) / 3\ J (s 



L (s) = In ( ( z r z + 1 + - (z r + z) + — — — - (z r + z + 

s s z } s V s / s 

(42) 



7 As we have already mentioned above it coincides with the action suggested in [Q, |) 
(see also 

8 Here we take (B — 0; N T = 2k) as an example, other cases can be easily obtained after 
the substitution of z — > — z or z r — > — z r . 
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Due to the fact that N% — > oo, the main contribution into £ comes from 
s ~ 0, so 

-M.^ln(^)-^ + 4 + 0( S 3) (43) 



with 



ft = + l) ; (44) 



and 

1 z r + z 



V = o7^- (45) 



3 z r z -\- . 

As far as N% — > oo, the function £ may be computed by the steepest 

descent method, with good accuracy. The saddle point so = so (x) is given 
by the equation 

dL(s ) 3 s / 2 x x /^x 

^r"3""^ +0(s » ) = i (46) 

so 

s ~ 4ft 0r? - X) (47) 

and we obtain 

- £ (s ) ~ -L (s ) + Sox ~ -3ft (x ~ ^Vj [X ~ \vj \ (48) 
Therefore we can finally write 

- 3 (x) = 3 (ft - A)) t + QvPoX + const (49) 
From the conditions 

- S ( X )' = 6 (p - ft) x + 677/3 ; -5 (*)" = 6 - ft) ; (50) 

we see that S (x) has the minimum for (3 < ft at 

*- = « ,s Tz^ft- (51) 



This solution (^T|) dominate for Xo < 4 , but at points xo — ±2 corresponding 
to 

p = /3 c = Po(l-l ; (52) 
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the regime changes. For (3 > f3 c the value Xmin keeps to 4, so 

Xmin = min { 4 ; xl) ; sign (xmin) = sign (77) (53) 



As it is seen from (|52"D, 'toy' fermions contribution leads to a shift of the 
effective couplingF]. For all m r , both fto and |f?| increase with increasing m. 
This leads to decreasing f3 c = Pc(y) at small m , but for m > T it starts 
to slowly increase to the asymptotic value (3 C (00) = |. 

An average value of (Q) of some Q = Q (cos |) can be computed simply 

as 

J 4 " Q (cos f) exp {s cos f } d/A / L(s) ( d 



^ ' J 4 "exp{scosf}ci/2 ) I \ds 



(54) 
with 

so - 4/3 ( -n - Xmm } (55) 



where Xmm is defined by (pal). 

In particular, taking into account that 

1 x 



2 r -COsf/ V 5s / S V 2 / (SfS+j)' 

(56) 

we obtain for (V^) 

^ = (z + (-l)^cosf) + ^ (z r + (-l) B cosf) (57) 
which for the case considered (|57|) gives for the even iV T 

<^> = (l - %T) (58) 



11 T 

and for odd N T 

9 It can be shown that such shift becomes a bit larger for a number of flavours n/ = 2, 
but leaves (3 C > 0, therefore does not wash out the phase transition. 



m + (1 _ ) (59) 
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It is easy to see that at critical point (where |x mm | — 2) both ( |5Tf ) and (5D) 
became comparatively small but still remain finite. 

A very important feature of the chiral condensate is the behavior of 
with decreasing quark mass m. Excellent fit of experimental data (in 
the symmetric phase of QED4) in the fermion mass range 0.01 — 0.06 is given 
in 



19 1 by the relation 



1 1 

(■4>4>) x ln e ~ m; e ~ - (60) 

(W 4 

The above roughly agrees with (]58|) that predicts that (4>ip) ~ 7p. It easy to 
check that the data in |2(J (obtained for 577(3) group with flavors number 
rif = 2 on lattice N T = 6 ) are consistent with Oipip) ~ 6am = ^ , which 
gives a reason to believe that fl58D indeed provides a reasonable description 
of the data. 



3 Discussion 

It is commonly assumed that, quenched lattice QCD studies using Wilson- 
Dirac fermions bring about large statistical errors in calculations involving 
very light quarks. Indeed, in quenched simulations, using Wilson fermions, 
lattice spacing effects were shown to be a major problem ||21|| . Fortunately, in 
the parameter area available for modern computers , sea quark contribution 
does not introduce substantial changes in MC data, except on occasions. For 
example, comparison of spectroscopy results obtained with dynamical and 
quenched fermions does not reveal any dramatic differences. Apparently at 
the parameter values of the simulation the sea quarks simply do not affect 
spectroscopy above the five to ten per cent level |S0|]. The measurements 
of valence observables in an almost- quenched system show the physics with 
light dynamical quarks to be very similar to quenched physics with the same 
p2| . The comparison between Wilson and clover actions showed that 



masses 



the results were agreeable |£J and no significant impact of unquenching was 



observed [24| on the hyperfine splitting as well. Therefore, the major part of 
MC data helps to trace the sea quark contribution, rather than obligatory 
makes us take the latter into account. 

Nevertheless, when the area of small quark masses and lattice spacing 
will be accessible in MC simulations, the situation presumably may change. 
Even with current limited computer capacity, there are some parameter areas 
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where the sea quark contribution is considerable. Indeed, in PSJ effects of 



dynamical fermions were singled out within the confinement phase, by com- 
parison of the quenched and full formulations of compact QED with Wilson 
fermions. It was established, that in the strong coupling limit (/3 = 0) the 
quenched theory was a good approximation of the full one but appeared to 
be in sharp contrast with it at /3 — 0.8. At such f3 the physics changes signif- 
icantly and the formation of metastable states was observed |25| presumably 
due to the presence of dynamical fermions. 

We hope that suggested toy model approximation will help not only to 
estimate the sea quark contribution at lattices with size and spacings un- 
accessible for MC experiment, but may help to qualify the parameter area 
where the presence of dynamical fermions leads to particularly appreciable 
effects. 



4 Conclusions 

Suggested approach to the problem of dynamical fermion effect estimation 
is based on the discard of all terms of l/£ order in the action. Therefore, 
only quark lines winding the lattice in the time direction are possible. The 
model allows an analytical solution for small values of lattice spacings and 
light quark masses on infinitely large lattices, where MC simulations eas- 
ily become prohibitively costly. For the vanishing chemical potential /j,, the 
fermion determinant remains real and can be incorporated into the mea- 
sure. Although such contribution is almost trivial, it significantly changes 
the phase structure of the model. It is believed that chiral symmetry vio- 
lated by Wilson term is restored in the continuum limit, however, no such 
restoration is revealed in the model at finite 2? and 

In fact, two issues, which are fundamental to understand the consequences 
of such approximation, have not been fully clarified: a remarkable difference 
between the results obtained on lattices with odd and even N T (see also p6| ) 
as well as very marginal alterity between periodic and antiperiodic border 
conditions. 

The first direct calculation of QCD properties for small quark masses was 



fulfilled long ago p7j . Substantial diversity between periodic and antiperi- 
odic border conditions was established. However, it was also pointed out 
that such difference might smear out with increasing lattice volume. As it 
is seen from (|4|) and (^6|) , alterity between periodic and antiperiodic bor- 
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der conditions formally does not disappear at any large N T and hardly will 
be still unobservable at higher orders in ^ and . As it follows from (|22|) 
(— 1) B factor that reflects the difference between periodic and antiperiodic 
conditions for even N may be absorbed into Z (iV)-transformation. There- 
fore, there is still a chance that even broken Z (2) -invariance may smear the 
mentioned difference, so an additional consideration of N = 3 is desirable, 

Another problem is that for large Wilson parameter r — > 1 (or m r — > 
oo ) the condensate < ipip > does not vanish at any finite temperatures. 
Moreover, at given approximation (in contrast to suggestion made in p5| ) 
chiral symmetry broken by the Wilson term can hardly be recovered by any 
fine-tuning of the bare parameters in the continuum limit. 



5 Appendix 

Expanding F n+1 



n+l 



I 


Qn+1 


b„ 


. 


. 


c n+l 


\ 




(-it 


Qn 


K-i ■ 


. 












c n-l 


q n -i ■ 


. 















o . 


• 12 


bi 




V 







o . 


■ Ci 


<h 


J 



(61) 



in elements of the first row 



d d 

Fn+l = Qn+lfn + b n ——F n+ i + C n+ i— ^n+lj 

ob n dc n+1 



(62) 



we obtain: 







n+l 



b n —F n+1 = -6 n c n / n _i- Y\ {~bj) 



(63) 



and 



where 



d 

C n+IT. F n+ i = —C n+ ib n+ if^ — Y\ (— Cj) 

n+1 j=l 



n+l 



fn — Qnfn-1 ~ ^n-l c n-l/n-2 — F n (b n — C n — 0) 



(64) 
(65) 
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and 

/f = j-j!- F n + l (66) 
dq n+1 dqi 

In a simple case of qj and the product bjCj independent of j, 

qj = g; 6 n = — ; c n = -A0 n (67) 

we can write 

/f = f, i (68) 

and 

- g/,- - A 2 /,-i = (69) 
Equation ( j69|) can be easily solved by introducing the generating function 

oo 

/(*) Y. f>' J 1 ( 7 °) 

i= n o 

This leads to 

f; Q W +1 - qfi** - aVjV" 1 ) = * 2 *) / (*) = ( 71 ) 

i=« ^ 

where n is an arbitrary fixed number. For example for n = 2 or n = 3 
from fl6~l~|) we easily find 

h = Q 2 + A 2 h = Q 3 + 2A 2 g (72) 
therefore, general solution may be written as 



f _ fn f z n 1 dz j 2A n cosh (n arcsinh^) ; n = 2k 
2ni J ± - g - X 2 z ~ \ 2A" sinh (n arcsinh ±) ) ; n = 2k + 1 



(73) 



Having collected everything we may finally write for F n 

cosh (n arcsinh ^) — ■^g— ; n = 2k; 
sinh (n arcsinh — e 2 ; n = 2A; + 1, 



(74) 



where 



e=n% (75) 
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